Rotational Sequences 


of Global Oscillations inside the Sun 


NASA 


Charles L. Wolff 

Goddard Space Flight Center, Greenbelt, MD 20771 


charles.wolffOgsfc.nasa.gov 


Rotational Sequences ot Global Oscillations inside the Sun 

by Charles Wolff 
POPULAR SUMMARY 


Many oscillation modes can exist inside the Sun and cover the entire globe. The most 
fundamental are modes trapped below the convective layers (i. e„ below the outer 30 /o ot 
a solar radius). These modes engage most of the Sun’s mass in organized motion and 
affect the deepest layers where its energy is generated. Beat periods between the rotation 
rates of the various modes can be many month to years. When the longer beat periods 
modulate solar luminosity, they become potentially interesting to the Earth climate 

problem. 

In this paper, I found clear evidence for the rotation rates of both gravity modes and 
Rossby modes. This proves they are excited in the Sun. There was also strong evidence 
for beat periods between the gravity modes. (Astronomers call these ‘‘g-modes and r- 
modes”. The Rossby waves of meteorology are just local approximations to the Rossby 

modes of a sphere.) 

All the low angular harmonic modes were detected--from #2 up to #7~by means of a 
very simple mathematical sequence which their rotation rates have to obey. Some 
additional modes up to #14 were also seen. Beyond that number, the observations do not 
have enough resolving power to reveal any higher modes. 

The observations used were the 54 year data set of daily measurements of thermal radio 
flux from the Sun. This rises and falls with other popular measures of solar variability 
such as the sunspot number and solar irradiance, so similar results would be expected it 
one of these data sets had been used. 

We learn two other things about the inside of the Sun in addition to detecting that the 
oscillations exist. The Sun’s average interior rotation period is measured to be 29.1 days, 
apparent from Earth. This is longer than the well known 27 day surface rotation period at 
low latitudes. Secondly, the 29.1 day period and theoretical properties of the modes 
require that the wave energy of the global modes is being reflected very deep in the Sun 
but not from the very center. The reflection layer is at 1 8% of the solar radius, which 
happens to be the same place where others have already measured a pronounced anomaly 
in sound speed. Something major is happening at 18% of the solar radius that 
theoreticians have yet to explain. 
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ABSTRACT 


A very simple mathematical sequence is detected in a half century of thermal 
radio dux from the Sun. Since the only known physical cause of the sequence is 
global oscillations trapped in the nonconvecting solar interior, g-modes and prob- 
ably r-modes are active. . If so, their rotation frequencies are detected and some 
previously reported difference frequencies are confirmed with high confidence. 


All angular harmonics for 2 ^ l ^ 


7 are detected as well as some others up to 


the limit £ <C 14 resolvable by the observations (a Fourier spectrum of the 10.7 


cm flux time series). The mean sidereal rotation of the nonconvecting interior is 
428.2 nHz as averaged by g-modes and 429.8 nHz by the r-modes, indicating that 
g-mode energy is a bit more centrally concentrated. Helioseismology measures 
such rotation rates near 0.36R (R = solar radius), so the global modes would have 
about half their kinetic energy above and below that level. This, and the known 
log(r) energy dependence of most modes implies that these oscillations are signif- 
icantly reflected near 0.18R, the same level at which sound speed measurements 
display a maximum departure from theoretical models. 


Subject headmgs: Sun: oscillations-Stars: oscillations-Sun: activity Sun 
interior 


1. INTRODUCTION 

Linear oscillation modes ?(r )e M trapped in stellar layers that rotate slowly at constant 
rate ft are well understood (Cox 1980; Unno, el ai 1989). In co-rotating coordinates, the 
oscillation frequency a and spatial variation, 


£ = rfyY + rfiVY - hr X VI', 


( 1 ) 
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depend on t.he radial and angular harmonic numbers (n, l , m), spherreal harmonic functions 
y = 4>), and the radial eigenfunctions /.ft )- Only one V/“ is needed to describe 

the dominant part of the motion which is consequently a running wave in longitude, . 
exp[i(m* + <rt)]. For r-modes the Coriolis force and h tern, dominate (Provost, Berthomieu, 
St Rocca 1981; Smeyers, Craeynest, k Martens 1981). For g-modes, buoyancy is the main 
force and the /, term is neglected (Cowling 1941). Christensen-Dalsgaard k Berthomieu 
(1991) reviewed g-mode properties and Wolff (2000) the r-modes. 


To first order, most g-modes and low radial harmonic r-modes have rotational 
properties that are particularly simple for standing modes (§2), the subject of this paper. 
This signature of their presence (eq. [2]) gives strict tests to apply to obseivations. Both g 
and r-modes would be trapped in the noncouvecting solar interior which rotates like a solid 
body with only a few percent deviation (Eff-Darwich, Korzennik, k Jimenez-Reyes 2002). 
This justifies using theory based on uniform rotation for early searches for the rotational 
signature. In §3, the signature is detected with high confidence in three different parts of 


the 10.7 cm flux spectrum. 


2. ROTATION OF GLOBAL MODES 

For r-modes, a point of constant phase drifts backward relative to a rotating stellar 
fluid. But relative to inertial space the mode turns forward at the rate (Saio 1982) 


where Us - Q/2tt is the fluid rotation rate in Hertz and the drift constant v d = 2 v s for 
r-modes. Since this is independent of m, the standing r-mode (sum of the +m and -in 
mode) rotates at the same rate. The inertial rotation of a standing g-mode is also given by 
equation (2) but now (Wolff 1974) u d = u.J/K from the classic solution of Ledoux (1951) 
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Where J = f dm(2hh + /|), K = / <M/«* + A'W + HI, rfm = *”*>*’ and " iS ^ 
density. Asymptotic g-modes (n ; » P) have J/K = 1 but, for other g-modes, J/K most 

be computed. 

A standard solar model was used (Guenther, et al. 1992) with its buoyancy frequency 
(Brunt- Vaisala) arbitrarily reduced for r < 0.21 R until reaching zero at r = 0.15R (R - 
solar radius). This simulates the main effect on these oscillations of an assumed partial 
mixing near 0.18R where sound speed measurements show a persistent discrepancy with 
interior models (Basil et al. 2000; Kosovichev et al. 1997). The steeply declining buoyancy 
frequency excludes oscillation energy from the inner core where most radiative damping 
would occur, making g-mode excitation more likely. Many have discussed mixing part or 
all of the core, e. g„ Richard & Vauclair (1997); Brun, Turck-Chieze, & Morel (1998). 

A sampling of g-mode solutions was computed and the resulting values of J/K are shown 
in table 1. Their departure from unity (the asymptotic case) closely follows the empirical 

formula 

J/K = 1 - 0.45(G/n) 3/i ( 3 ) 

where t = [t{i + 1)] 1/2 . Oscillations trapped between 0.18R and the convective envelope, 
0.713R, will have equal kinetic energy above and below the logarithmic midpoint, since 
energy is distributed closely proportional to log(r) for most g-modes in the Sun (Wolff 
1979) and also for r- modes. At this midpoint, 0.358R, the Sun’s sidereal rotation is 429 ± 2 
nHz from the curves given by Eff-Darwich, Korzennik, & Jimenez-Reyes (2002). The rate 
has a slight positive slope with radius. Since the effective rotation sensed by an oscillation 
is the energy-weighted mean, v s in equation (2) should be close to 429 nHz. 
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3. ROTATIONAL SIGNATURES IN SOLAR ACTIVITY 

The 10.7 cm Radio Flux Index F l0 is available at www.drao.nrc.ca/icarus/www from 
Dominion Observatory. The ~ 20,000 daily values from 1947 Mar 5 through 2001 Nov 18 
were compressed to about 4000 by taking 5 day averages. This widely used index vanes 
with the 11 year cycle and the monthly rotation. Any true frequency v t influencing the 
index appears in the Fourier spectrum with aliases u t ± w m where i is a small integer and 
modulating frequencies «, m are (11 yr)" 1 , (1 month)- 1 , and any other strong modulator. 
This clutters the spectrum with many lines that look like noise. To reduce long period 
aliases, the deviation from the running mean F - a one year boxcar average-was used. 

Figure la shows the Fourier spectrum of (F l0 - F)F~*. The decline below 30 nHz is 
not solar but is due to subtracting F. The feature between 400 and 470 nHz is sometimes 
loosely attributed to active regions rotating across the disk but the fastest such rate 
apparent from Earth is 435 nHz. At least part of this bulge is unexplained. The rise in 
amplitude below 60 nHz may be due to overlapping sequences of beats between g-mode 
rotation rates mostly following equation (2) (Wolff (1983), hereafter ’’paper I”). The one 
beat sequence in that paper that overlaps no other is strongly confirmed on Figure lb by 
the relevant part of the F l0 spectrum. Vertical lines mark the differences between the l 
= 2 case and the other frequencies in table 1 of paper I, extended for l > 11 using the 
same equation. These theoretical lines are rotational beat frequencies. They mostly hit 
tall lines in the F 10 spectrum and every line for 8 is detected. Anyone assuming that 
the observed spectrum is mostly noise bears the heavy burden of explaining why the same 
theoretical sequence is detected in both F 10 (Fig. lb) and in sunspots (paper I) when only 
25% of the latter data is contemporaneous with the former . Noise doesn’t behave that way. 


The well known binomial expression 
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P(n, k) 


!>': 


7^(1 -P) 


n—k. 


( 4 ) 


k\(n — k)\ 

gives the probability P that k frequencies out of a total group of n will accidentally agree 
with observation if p is the probability that any random frequency agrees. To use this we 
need a success criterion for which p is a constant. The dashed curve, 1.2M, on every figure 
is based on M — a 20 nHz wide running mean of the spectrum. It turns out that 32.4% of 
the spectral points lie above this curve in any broad segment of the spectrum. By taking 
the dashed curve as the threshold for detection, one can use p = 0.324 everywhere. Of the 
12 resolvable lines on Figure lb, nine hit the spectrum above this curve making P(12, 9) = 
0.0027, which is an impressive 3a detection of the beat sequence. 


When applying equation (2) to g-modes, = 381 nHz as determined in paper I which 
stated it more generally as (762 nHz)/S where S (= 1 or 2) is the rotational symmetry of 
the modes. In all subsequent work, S = 2 has been used, meaning that a family of g-modes 
for a given l has two hot longitudes, 180" apart at low latitudes (Wolff 1974). Thus, 

2 („, - are the expected beat frequencies with the l = 2 g-mode, including the small 
deviations measured in paper I for l i 3. Finally, note that u. cancels from the difference 
frequencies so paper I actually measured u d and not the stellar rotation rate. 

In a sufficiently long data set, the rate at which a hot g-mode longitude faces Earth, 
2(y t - 31.69 nHz), should also appear. These lines should be weaker than the beats since 
detectible change at the surface is likely only after some area at the base of the convection 
zone has been heated enough to stimulate an overturning. Beats do this better since two or 
more modes are contributing at the same location. Searching close to the expected 429 nHz 
(see end of §2), I found that equation (2) with u s = 428.2 nHz fits the spectrum. Figures 
1c and Id mark these g-mode rates. Of the 13 predicted lines resolvable by the data, 8 
hit the spectrum above the dashed curve. Since a free parameter {u a ) always allows one 


observed line to be matched, I discarded one successful line from the sample obtaining 
P ( 12j 7) = 0.042. This 2a confidence level lends support to the above 3a detection of the 
same sequence by the beat frequencies. The overall probability that agreement on figure 1 
is accidental is P(12+12, 9+7) = 5 x 10 4 . 

The r-modes assist convection not with heat but by supplying vorticity and a zonal 
velocity gradient. This mechanism, reinforced by the £ = 1 mode, favors the synodic 
frequencies, u t - 31.69 nHz, which are shown on figure 2 using i/ a = 429.8 nHz in equation 
(2). Nine of the 13 lines are detected including all for £ ^ 7 . Again discarding one 
successful line gives an accidental probability, P(12, 8) = 0.013, which is a 2 \o detection 
even without correcting for the interference at £ = 9 and possibly £ = 8 from the dotted 
triplet. (Closely spaced lines can either depress the spectrum if out of phase or enhance it if 
in phase.) The triplets are centered on the sidereal rates u t . Evidence that these are aliased 
into at least triplets as they beat with the £ = 1 modes will be discussed elsewhere. Wolff & 
Hoegy (1989) suggested 1.5 nHz for such a rate and that is the splitting drawn on Figure 2. 

4. SUMMARY AND IMPLICATIONS 

Few realize how hard it is to fit data with a large number of rates from equation (2) 
having only one free parameter, i/,. (v d was already known from theory for r-modes and 
from earlier work for g-modes.) The g-mode sequence was detected with especially high 
confidence: accidental probability, 5 x 1(T 4 . Unfortunately, seeing high confidence levels 
has little impact existing opinions. But it must strain belief that we have been uncommonly 
lucky four times (!) with the same sequence (figs, lb, lc-d, 2, and paper I). As the figures 
showed, lines expected for £ < 8 were often prominent in the F 10 spectrum and exceeded 
the detection threshold 15 out of 17 times. Lines for higher l were detected less often (11 
of 21) but, still more than the 32% random expectation. The weakness for higher £ is not 
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surprising since less information about smaller horizontal scales can make it up through the 
convection zone to affect the surface. For example, oscillation energy declines roughly as 
(0.71R/r) 2 ' for the higher t modes. Evidence for l = 1 modes would require more space 

than is available here. 

A standard solar model was used with the assumption that oscillations suffer significant 
reflection near 0.18R where the sound speed is anomalous. The assumption found support 
when values of i/, near the expected 429 nHz gave good fits to the spectrum. Also, reflection 
of energy near 0.18R roughly doubles g-mode oscillation periods, partly explaining why 
their expected periods have not been detected (Palle 1991) even though their rotation 
rates are seen. Another implication follows from the ratio v d /v a = 381/428.2 determined m 
§3. This means that J/K = 0.890 and, by equation (3), the typical g-mode has a radial 
harmonic number n 2 2.W which is nearly in the asymptotic range. Finally, only the 
lowest degree g-modes (( ^ 3) deviate from equation (2) (table 1, paper I). No r-mode 
deviates except perhaps l = 1. Thus, second order terms neglected in deriving equation (2) 
are too small to be sensed in a half century of the 10.7 cm flux index except for the very 


lowest angular states. 
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Table 1. CORIOLIS EFFECT 
ON <7- MODES 


l n J/K i n J/K 


3 5 .7142 

3 10 .9090 

3 15 .9526 

3 20 .9699 

5 5 -4556 

5 6 .5738 

5 8 .7243 

5 10 .8087 


5 

20 

.9364 

5 

25 

.9549 

10 

10 

.5236 

10 

15 

.7430 

10 

20 

.8372 

10 

25 

.8803 

15 

15 

.5576 

15 

20 

.7109 

15 

25 

.7924 


5 15 .9020 



Frequency (nHz) 


Fig. 1 -- The spectrum of a half century of the 10.7 cm solar radio flux index and expanded segments, (a) The mound peaking at 430 nHz is 
associated with solar rotation rates and the rise below 60 nHz may be due to g-mode beats, (b) The sequence of rotational beats between the 
L = 2 g-mode and modes of higher L is strongly detected in the spectrum up to L = 8 and probably beyond. The detection threshold 
(dashed curve) lies above 2/3 of the spectrum, (c & d) The corresponding rates at which the g-mode hot longitudes face Earth are also 
detected, but less firmly because they are weaker, as expected. The 10.7 cm record is too short to resolve any mode for L > 14, indicated by 
the flat rectangle ending at L = infinity. 









Fig. 2- Part of the same observed spectrum on Figure 1 . Theoretical lines are now the synodic rotation rates of r-modes. Again, all low 
angular harmonics are detected as well as some others up to the resolution limit. The L = 9 case suffers interference from one of the dotted 
triplets located where r-mode beats (and aliases) with the L = 1 modes should be located. 





